Constraints on inflation from CMB and Lyman-a forest 



O 
O 
(N 

t^: 
< 

(N 
> 

o 
cn 
o 

o 

:^ 

o 
U 



X 



S. Hannestad ^ 

NORDITA, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
S. H. Hansen ^ 
NAPL, University of Oxford, Kehle road, 0X1 3RH, Oxford, UK 
F. L. Villante ^ 

Dipartimento di Fisica and Sezione INFN di Ferrara, Via del Paradiso 12, 44^00 Ferrara, Italy 

A. J. S. Hamilton ^ 
JILA and Dept. Astrophys. & Planet. ScL, Box 440, U. Colorado, Boulder CO 80309, USA 

Abstract 

We constrain the spectrum of primordial curvature perturbations by using recent 
Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. Specifi- 
cally, we consider CMB data from the COBE, Boomerang and Maxima experiments, the 
real space galaxy power spectrum from the IRAS PSCz survey, and the linear matter 
power spectrum inferred from Ly— a forest spectra, where we for simplicity assume the 
absence of appreciable covariances. We study the case of single field slow roll inflation- 
ary models, and we extract bounds on the scalar spectral index, n, the tensor to scalar 
ratio, r, and the running of the scalar spectral index, di^k, for various combinations 
of the observational data. We find that CMB data, when combined with data from 
Lyman— a forest, place strong constraints on the inflationary parameters. Specifically, 
we obtain n « 0.9, r < 0.3 and dy^k ~ 0, indicating that hig n, big r models (often 
referred to as hybrid models) are ruled out. 

PACS: 98.62.Ra, 98.65.-r, 98.70.Vc, 98.80.Cq 

1 Introduction 

Inflation is generally believed to provide the initial conditions for the evolution of large 
scale structure (LSS) and the cosmic microvi^ave background radiation (CMB). The garden 
of inflation offers a bounty of models, each of which predicts a certain power spectrum of 
primordial curvature perturbations, V{k), a function of the wavenumber k. This power 
spectrum can be Taylor-expanded about some wavenumber fco and truncated after a few 
terms [1] 
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in which the first term is a normalization constant, the second is a power-law approximation, 
with the case n = 1 corresponding to a scale invariant Harrison-Zel'dovich spectrum, and 
the third term is the running of the spectral index. 

An important class of models is given by single field slow roll (SR) inflationary models, 
which can be treated perturbatively. The properties of SR models are well known (see e.g. 
ref. [2] for review and a list of references), and we will here classify the different models 
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by the 3 parameters (n,r, 9infe)j where n = dlnV/dlnk\k=ko + 1 is the scalar spectral index 
at the pivot scale k^, the parameter r is the tensor to scalar perturbation ratio at the 
quadrupole scale, and di^k = dn/dlnk\k=ii.Q. The reason for using these 3 variables instead 
of the normal SR parameters (e,r/, (defined in the appendix) is simply that the former 
3 variables are more closely related to what is measured from observations. 

In SR models the tensor spectral index and its derivative can be expressed [2, 3] 
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The factor k in the above equations depends on the model, and in particular is different 
for different [4]- In our analysis we use the parametrisation from ref. [5], which for the 
models considered here means 5 < k < 7. 

The different SR models are traditionally [6] categorised into 3 main groups according 
to the relationship between the first and second derivative of the inflaton potential, and 
they are distributed in (n, r) and (n, 9infc) space as in Fig. 1 (see details in the appendix 
or in ref. [7]), where the dashed lines on the borders between the different models are the 
two attractors found in ref. [8], r = and r = —K,{n — 1). For large values of n — 1 and 
r, SR models "naturally" predict sizeable deviations from a power law approximation (i.e. 
d\ak 7^ 0), so when comparing to observations it is important to include the third parameter 
d\ak^ in addition to the parameters n and r commonly considered. 




0^2 



-0.1 



0.60 0.70 0.80 



90 1 .00 1.10 1 .20 
n 




0.0 



0.60 0.70 



30 0.90 1.00 1.10 
n 



Figure 1: The various slow-roll models in (n, r) and {n^diak) space. The dashed lines are 
the two attractors (here we have used k = 5 which is a typical value for a flat universe with 
a large cosmological constant). The figure (ra, c^infc) is for = 0. The hatched regions move 
up and down by inclusion of the third derivative, 7^ (see Eq. 10 in the appendix). 

The purpose of this paper is to ascertain what constraints currently available observa- 
tional data can place on the parameters {n,d\-ak^f) of SR inflation. Of course the universe 
could lie outside the hatched regions of Fig. 1, indicating that one should look beyond SR 
inflation; but it seems reasonable to explore the simplest models in the first instance. In 
ref. [7], CMB data were used to put bounds on inflationary parameters, and indications 
were found towards a negative bend of the primordial power spectrum (9infc ^ 0). However, 
the scales probed by CMB are limited, and the constraint on d\-ak was weak. 

In the present paper we extend the analysis of [7] by including information not only from 
the CMB, but also from the linear galaxy power spectrum measured from the IRAS Point 
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Source Catalogue Redshift (PSCz) survey [9], and from the linear matter power spectrum 
inferred from the Ly— a forest in quasar spectra [10]. Even though it is not completely clear 
that the correlations in the 3d matter power spectrum are negligible, we treat them as such 
for simplicity. The size of and effect from such covariances might be essential, and should 
be considered more carefully in a future investigation. As we will show, these data, probing 
different ranges of scales, yield much tighter constraints on the inflationary parameters than 
CMB data alone. 



Recently observational data on both the CMB and LSS have improved dramatically. The 
balloon-borne experiments Boomerang [15] and MAXIMA [16] have determined the CMB 
angular power spectrum beyond the first acoustic peak, whose position near I ~ 200 in- 
dicates a fiat universe, hence seemingly confirming the inflationary paradigm. The two 
experiments, together with COBE-DMR [18], provide a high signal-to-noise determination 
of the CMB power spectrum for 2 < / < 800, thus testing the structure formation paradigm 
on scales roughly of the order k < 0.1 h Mpc~^. 

LSS data, probing smaller scales, provide complementary information. In this paper we 
consider the real space galaxy power spectrum from IRAS Point Source Catalogue Redshift 
(PSCz) survey [19] which probes scales in the range k = 0.01 — 300 h Mpc~^ [9, 17]. In 
order to avoid problems with the interpretation of non-linear effects, we use data only at 
scales k <0.3 h Mpc~^ 

Further, we consider the information on the linear matter spectrum at redshift z ~ 3 
which can be inferred from the Ly-a forest in quasar spectra. The fact that the nonlinear 
scale is smaller at higher redshift makes it possible to probe the linear power spectrum to 
smaller scales {k ~ 0.2 — 5 h Mpc~^) than are accessible to galaxy surveys at low redshift. 
In this paper we use a recent determination [10] of the linear matter power spectrum from 
the Lyman-a forest at redshift z = 2.72, based on a large sample of Keck HIRES and Keck 
LRIS quasar spectra. To convert the measured power spectrum of Ly-a flux into the matter 
power spectrum, [10] apply a correction factor obtained from A'^-body computer simulations. 
To allow for possible systematic uncertainty [20] in this correction factor, we will repeat the 
analysis excluding the data at the smallest scales. 

2.1 Data analysis 

In order to investigate how the CMB, PSCz and Ly-a data constrain the SR parameter 
space (n, r, we performed a likelihood analysis of the data sets from COBE [18], 

Boomerang [15] and MAXIMA [16], together with the decorrelated linear power spectrum 
of PSCz galaxies for k<O.S h Mpc'^ [9, 17], and the Ly-a data from Table 4 of ref. [10]. 
The likelihood function is 



2 The data 
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For the CMB data 
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while for PSCz and Ly-a data 
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the sum being taken over pubhshed values of band-powers Ci^i and Pfe^j. The quantity 9 is 
a vector of cosmological parameters, taken here to be 



The parameters are: the matter density Q^; the baryon density the Hubble parameter 
Hq; the optical depth r to reionization; the normalizations QcMB, QpsCz, and Qhy-a of 
the CMB, LSS, and Ly-ct power spectra; and the inflationary parameters n, r, Sin^. We 
have assumed that the universe is flat, JIa = 1 — ^m, as predicted by standard inflationary 
models. For all the figures we marginalize over all other parameters. We consider only two 
values for the baryon density, Cli,h^ = 0.019 and 0.030, however, as we will see, the results 
are very similar in those two cases, suggesting that the results will be similar if allowing 
Q^h^ as a free parameter. 

The CMB, LSS, and Ly-a data are all subject to uncertainties in their overall normal- 
izations. The CMB groups quote estimated calibration errors for their experiments, which 
we account for by allowing the data points to shift up or down, by 10% for Boomerang [15], 
and by 4% for MAXIMA [16]. Because of uncertainty in the linear galaxy-to-mass bias for 
PSCz galaxies [17], we conservatively treat the normalization Qpscz as an unconstrained 
parameter. Ref. [10] quotes uncertainties in the overall normalization of the matter power 
spectrum inferred from the Ly-a forest, but these uncertainties are based on simulations 
with 5infc = 0, so to avoid possible bias we again conservatively treat the normalization 
Qhya as an unconstrained parameter. 

For simplicity all data points have been treated as uncorrelated in the likelihood func- 
tions, eqs. (5,6). For the CMB data, correlations between estimates C;^j of angular power 
at different harmonics I are induced by finite sky coverage, but in practice the CMB teams 
quote band-powers at suflicicntly well-separated bands of I that the correlations are prob- 
ably small. For the PSCz data, the published band-powers are explicitly decorrelated. For 
the Ly-a data, the covariances between estimates of the flux power spectrum are small, 
according to Fig. 12 of [10], and this may translate into small statistical covariances in the 
inferred matter power spectrum. As mentioned earlier, it is not completely clear how good 
this translation from flux power to matter power is, and we leave this question for future 
investigation. 

We have chosen the pivot scale in Eq. (1) as kQ = 0.05 Mpc~"^. This choice is made for 
convenience, since ko = 0.05 Mpc^^ is the scale at which wave-numbers are normalised in 
the CMBFAST code. Our results are independent of the value of ko. 
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Figure 2: The 1 and 2 a allowed regions for the two slow roll parameters n and di^k- The left 
panels assume a BBN prior on Q^h'^ = 0.019, whereas the right panels are for Qbh"^ = 0.030, 
the value which best fits the CMB data. The top row is for CMB data alone, the middle 
row is for Lyman-a data alone, and the bottom row is for the combined analysis. 
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2.2 Results 



The analysis of the constraints from CMB data was presented in detail in ref. [7] (see 
also [21]). Here we extend the analysis in [7] by including the reionization optical depth, r, 
as a free parameter. This is potentially important since there is a well-known degeneracy 
between r and the scalar spectral index n. However, the analysis turns out to prefer models 
with T « 0, so including r leaves the main conclusions of [7] essentially unchanged, as 
can be seen in the upper panels of Figs. 2 and 3. We recall here the main conclusions 
from [7]: (1) we allow the primordial power spectrum to bend, dij^k 7^ 0, then CMB data 
do not constrain the tensor to scalar perturbation ratio r; (2) if we assume a BEN prior, 
Q^h^ = 0.019 [24], then CMB data favour a negative bend, di^k < 0) corresponding to a 
bump-like feature centered at scales k ^ 0.004 Mpc~^. 

Let us now discuss the information on the inflationary parameters which emerges from 
PSCz and Ly-a data. The first observation is that PSCz data do not provide relevant 
constraints in the (n, Oin^) plane. This is easily understood, because the PSCz data at 
large scales (say k < 0.03/tMpc~^) have large errors and therefore play little role in the 
evaluation. This implies that PSCz data effectively span only one decade in k. Considering 
that the overall normalization of the data is taken as a free parameter, it is evident that 
one cannot obtain strong constraints from such a small range of scales. 

The situation is quite different with Ly-a data, which have small error bars and span 
almost two decades in k. This is shown in Fig. 2, where we present the allowed regions 
corresponding to 1 and 2 a (we define the la region as 

< Xmin + 2.3 and the 2a region 
as < Xmin + 6-17) for both Qhh'^ = 0.019 as suggested by BBN (left column), and 
Qfe/i^ = 0.03 as suggested by CMB (right column). The top graphs are from CMB data 
alone, the middle graphs are obtained using Ly-a data alone, and the bottom graphs are 
from the combined analysis. 

It is straightforward to understand how Ly-a data select the allowed regions shown in 
the middle panels of Fig. 2. At the small scales probed by Ly-a data, the theoretical linear 
matter power spectrum P{k) is roughly proportional to 

P(/e) ocln2(aA;/r/i)/e"'-^ (8) 

where a = 2.205, F = i^mh and k is expressed in Mpc~^ [25]. The parameter n' = 
dlnVika) / din k + 1 is the effective spectral index of primordial density perturbation at the 
scale ka ~ 3.7(0^/1^/0.4)^/^ Mpc~^ representative of Ly-a data ^. On the other hand the 
Ly-a data, as discussed in [10], are well fitted by a power law, P{k) ~ k", with spectral 
index = -2.47 ± 0.06. This means that 

^hpM = (n' - 4) + ^ , ,^ = -2.47 ± 0.06 . (9) 
dink ^ ' ln{aka.lTK) ^ ' 

If we consider that ^\)k? < Th < 1, the previous expression can be directy translated into 
a bound for the effective spectral index of primordial density perturbations. We obtain 
n' ~ 0.75 — 1.1, which roughly corresponds to what is shown in the middle panels of 
Fig. 2 for Oinfe = 0. It is also easy to understand the observed correlation between n and 
^infc. The effective spectral index n' can be expressed, as a function of n and di^ki by 
n' = n + ln{ka/kQ)d\^k- It is thus evident that the allowed region in the plane (n, 9infc) lies 
roughly along lines of constant n'. 

®The observational units for wavenumbers are (kms~^)~^. The conversion to Mpc~^ is model-dependent 
since it requires the evaluation of the Hubble constant H{z) at redshift z = 2.72. 
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It is important to note that the regions constrained by Ly-a data are "orthogonal" to the 
regions constrained by CMB data. This means that the combined analysis (CMB+Ly-a) 
gives much stronger bounds than any of the two alone. Moreover, since CMB data provide 
a bound on Q.mh'^i the degeneracy between n' and Th in Eq. (9) is removed, and therefore 
the Ly-a observation of P{k) oc /j-2.47±o.06 (directly translated into a constraint on 

n' . This is clearly shown in the lower panels of Fig. 2, from which one obtains the following 
conclusions: (1) the combined analysis strongly indicates that the bend of the primordial 
power spectrum is close to zero, being —0.05 < c^infe < 0.02 at the 2a level; (2) the spectral 
index n is constrained to be fairly close to n = 0.9. One notes that the limits obtained do 
not crucially depend on the assumed values of $7;,, indicating that the inclusion of Lyman-a 
data avoids the problem with CMB data alone, that different fi^ extend the allowed range 
of the spectral index beyond n ~ 1 [21]. However, one should note that, due essentially 
to CMB data, the goodness of the fit is quite sensitive to the assumed value of 0;,, the 
being substantially smaller for high 17;,. 
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Figure 3: The 1 and 2 a allowed regions in the (n, r) plane. The left panels assume a BBN 
prior on i^^h'^ = 0.019, whereas the right panels are for Q^/i^ = 0.030, the value which best 
fits the CMB data. The top panels are for CMB data alone, while the bottom panels are 
for the combined analysis. Hybrid models are to the right of the full line. 

In Fig. 3 we present the constraints obtained for the remaining SR parameter r. Specif- 
ically, we show the 1 and 2 a allowed regions in the (n, r) plane, for both low and high 0^. 
The top panels are from CMB data alone, while the bottom ones are from the combined 
analysis. It is clear that in the combined analysis r is constrained to be smaller than about 
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0.3 at 2 a. This is substantially different from the analysis of CMB data alone, where no 
constraints on r could be obtained, since a strong bend, e.g. c?infc = —0.2, could allow the 
tensor component to be big [7]. It is worth noting, that Ly-a data don't probe r directly; 
they fix n to be less than 1, and di^^k close to zero. In this way strong constraint on r 
can be obtained from CMB data. This result is extremely important when we compare to 
theoretical models, since different classes of models predict different relationships between 
the scalar spectral index n and the tensor to scalar ratio r (see Fig. 1). The combined 
analysis, indicating n ~ 0.9 and r < 0.3, favours small field models and seems to exclude 
hybrid models. This is clear from the lower panels of Fig. 3, where the big n, big r (often 
referred to as hybrid models, to the right of the full line), are excluded at 2a. ^ 

An important caveat to this analysis concerns possible systematic uncertainties in the 
inference of a linear matter power spectrum from the Ly-a data. What [10] measure directly 
from observations is the power spectrum of transmitted Ly-cc flux. The conversion to a linear 
matter power spectrum involves a fairly large, scale-dependent correction which [10] extract 
from collisionless computer simulations, with the Ly-a optical depth taken proportional to 
a certain power of the dark matter density [13]. It has been suggested [20] that the relation 
between baryonic and dark matter densities could introduce significant uncertainty in the 
flux-to-mass correction at small scales. Specifically, pressure effects cause the baryonic 
density 5^ to be smoothed compared to the dark matter density 6^, an effect that can be 
parametrized as (5f,(k) = exp[— (A;//c/)^](5d(k), with comoving filter scale kf ~ S5hMpc~^ 
[14]. The procedure considered by [20] is to treat the filter scale fcj as a free parameter 
constrained only by the shape of the power spectrum of Ly-a flux. However, [10] argue that 
treating kf as a free parameter is overly pessimistic, and that li kf takes values suggested 
by hydrodynamic simulations, then the effect on the power spectrum is minor. 

To allow for the possibility that the systematic errors are underestimated at small scales, 
we repeated the analysis neglecting the last 3 data-points from the Ly-a data [10], both for 
5infc = and for 9infc free. The upper panels in Fig. 4 show the case where di-^k = 0- Here, 
the results are essentially unchanged by the removal of the data points. However, the lower 
panels show the full case. Here, the tight constraint on r disappears completely, and the 
constraint on n is significantly weakened. The reason is that the small scale data points 
are the most important for constraining di^k- When these points are removed, dink is not 
nearly as tightly constrained as before, and as seen for the case where only CMBR data is 
used, a large r can be compensated by a negative bend of the spectrum. 

Thus, the very tight constraint derived above depends on the correctness of the small 
scale Ly-a data. Therefore, it is highly desirable that a better understanding of the possible 
systematic errors in determining matter power spectra from Ly-a forest observations is 
developed. 

Another worry is that the error ellipses for CMB and Ly-a are only marginally over- 
lapping. This might indicate that the two data sets are mutually inconsistent, and that 
other physical effects should be taken into account. For the Ly-a data one might think of 
massive neutrinos or warm dark matter, where both would suppress power on small scales, 
potentially allowing bigger n or bigger positive di^k- 

®We have used k = 5 (see Eq. 2), and for larger k the hybrid models move to bigger r. It is worth 
pointing out, that the classification of hybrid models [11] as in Fig. 1, and as the models to the right of 
the full line in Fig. 3 is oversimplified. In this paper we follow ref. [6], and by "hybrid models" we refer to 

potentials for which e < rj (see appendix for details). When considering more general potentials, or F-term 
hybrid inflation [12], more complicated behaviour results, and a simplified classification is impossible. 
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Figure 4: The 1 and 2 cr allowed regions in the (n, r) plane. These results are obtained from 
the combined analysis, neglecting the last 3 data-points from [10]. The left panels assume a 
BBN prior on Qf,h'^ = 0.019, whereas the right panels are for flf,h'^ = 0.030, the value which 
best fits the CMB data. The top panels assume di^k = 0, whereas the lower panels have 
dink as a free parameter. Hybrid models are to the right of the full line. 

3 Conclusions 

We have considered data from both CMB and LSS to place constraints on the parameters of 
single field slow roll inflationary models. We have found that, by combining CMB data with 
the power spectrum inferred from Lyman— a forest, one obtains strong constraints on the 
inflationary parameters. We obtain 0.8 < n < 1.0, r < 0.3 and —0.05 < di^k ^ 0.02 at 2a 
level. In the language of SR this means e < 0.03 and rj < 0.06, with the best fit model being 
small fields. This still leaves a large part of the SR parameter space open, but seems to 
exclude big n, big r models, often referred to as hybrid models. These constraints are much 
stronger than those from CMB data alone, and arise because the error-ellipse from Ly-a 
data is almost perpendicular to the error-ellipse from CMB data. Let us repeat, that we 
in this analysis for simplicity have assumed the absence of appreciable covariances between 
the reconstructed Ly-a mass powers. 
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A Notation 

We use the notation: r] = ae, where 

see [7, 21] for details. The notation with r] = M'^V" /V , used e.g. in [2], simply corresponds 
to the substitution a ^ a + 1. The 3 classes of SR models are small fields {a < —1), large 
fields (— 1 < a < 1), and hybrid models (1 < a). One finds 

r = 2Ke, ainfe = -2^2 + 8e^(2a-l) and n - 1 = 2e(a - 2) . (10) 
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Abstract 

We constrain the spectrum of primordial curvature perturbations by using recent 
Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. Specifi- 
cally, wc consider CMB data from the COBE, Boomerang and Maxima experiments, the 
real space galaxy power spectrum from the IRAS PSCz survey, and the linear matter 
power spectrum inferred from Ly— a forest spectra, where we for simplicity assume the 
absence of appreciable covariances. We study the case of single field slow roll inflation- 
ary models, and we extract bounds on the scalar spectral index, n, the tensor to scalar 
ratio, r, and the running of the scalar spectral index, d\nk, for various combinations 
of the observational data. We find that CMB data, when combined with data from 
Lyman— a forest, place strong constraints on the inflationary parameters. Specifically, 
we obtain n « 0.9, r < 0.3 and d\nk ~ 0, indicating that hig n, big r models (often 
referred to as hybrid models) are ruled out. 

PACS: 98.62.Ra, 98.65.-r, 98.70.Vc, 98.80.Cq 



1 Introduction 



Inflation is generally believed to provide the initial conditions for the evolution of large 
scale structure (LSS) and the cosmic microwave background radiation (CMB). The garden 
of inflation offers a bounty of models, each of which predicts a certain power spectrum of 
primordial curvature perturbations, V{k), a function of the wavenumber k. This power 
spectrum can be Taylor-expanded about some wavenumber ko and truncated after a few 
terms [1] 



lnV{k) = lnV{ko) + {n 



_^ /i- 1 dn 
^ ^ 2d\nk 



ko ko 



(1) 



in which the first term is a normalization constant, the second is a power-law approximation, 
with the case n = 1 corresponding to a scale invariant Harrison- Zel'dovich spectrum, and 
the third term is the running of the spectral index. 

An important class of models is given by single fleld slow roll (SR) inflationary models, 
which can be treated perturbatively. The properties of SR models are well known (see e.g. 
ref. [2] for review and a list of references), and we will here classify the different models 



^e-mail: steen@nordita.dk 
^e-mail: hansen@astro.ox.ac.uk 
^e-mail: villante@fe.infn.it 
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by the 3 parameters {n,r,dink), where n = dlnV /dlnk\k=kf^ + 1 is the scalar spectral index 
at the pivot scale ko, the parameter r is the tensor to scalar perturbation ratio at the 
quadrupole scale, and dink = dn/d\nk\k=ko- The reason for using these 3 variables instead 
of the normal SR parameters (e, ?7,^^) (defined in the appendix) is simply that the former 
3 variables are more closely related to what is measured from observations. 

In SR models the tensor spectral index and its derivative can be expressed [2, 3] 



r driT r 

riT = and — — - = - 

K d ink K 



1) + 



(2) 



The factor k in the above equations depends on the model, and in particular is different 
for different f^A [4]- In our analysis we use the parametrisation from ref. [5], which for the 
models considered here means 5 < k < 7. 

The different SR models are traditionally [6] categorised into 3 main groups according 
to the relationship between the first and second derivative of the inflaton potential, and 
they are distributed in (n, r) and {n,dink) space as in Fig. 1 (see details in the appendix 
or in ref. [7]), where the dashed lines on the borders between the different models are the 
two attractors found in ref. [8], r = and r = — K(n — 1). For large values of n — 1 and 
r, SR models "naturally" predict sizeable deviations from a power law approximation (i.e. 
dink 7^ 0)5 so when comparing to observations it is important to include the third parameter 
dink, in addition to the parameters n and r commonly considered. 




0.60 0.70 0.80 0.90 1.00 1.10 1.20 
n 




-0.1 



0.60 0.70 0.80 0.90 1.00 1.10 
n 



Figure 1: The various slow-roll models ill (77., r) ciiid. (?7-, ^njt 

) space. The dashed lines are 

the two attractors (here we have used k = 5 which is a typical value for a fiat universe with 
a large cosmological constant). The figure (n, dink) is for = 0. The hatched regions move 
up and down by inclusion of the third derivative, / (see Eq. 10 in the appendix). 

The purpose of this paper is to ascertain what constraints currently available observa- 
tional data can place on the parameters (n, dink^i") of SR infiation. Of course the universe 
could lie outside the hatched regions of Fig. 1, indicating that one should look beyond SR 
infiation; but it seems reasonable to explore the simplest models in the first instance. In 
ref. [7], CMB data were used to put bounds on infiationary parameters, and indications 
were found towards a negative bend of the primordial power spectrum (dink ^ 0). However, 
the scales probed by CMB are limited, and the constraint on dink 

was weak. 

In the present paper we extend the analysis of [7] by including information not only from 
the CMB, but also from the linear galaxy power spectrum measured from the IRAS Point 
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Source Catalogue Redshift (PSCz) survey [9], and from the linear matter power spectrum 
inferred from the Ly— a forest in quasar spectra [10]. Even though it is not completely clear 
that the correlations in the 3d matter power spectrum are negligible, we treat them as such 
for simplicity. The size of and effect from such covariances might be essential, and should 
be considered more carefully in a future investigation. As we will show, these data, probing 
different ranges of scales, yield much tighter constraints on the inflationary parameters than 
CMB data alone. 



2 The data 

Recently observational data on both the CMB and LSS have improved dramatically. The 
balloon-borne experiments Boomerang [15] and MAXIMA [16] have determined the CMB 
angular power spectrum beyond the first acoustic peak, whose position near I ~ 200 in- 
dicates a flat universe, hence seemingly confirming the infiationary paradigm. The two 
experiments, together with COBE-DMR [18], provide a high signal-to- noise determination 
of the CMB power spectrum for 2 < / < 800, thus testing the structure formation paradigm 
on scales roughly of the order A-' < 0.1 h Mpc~^. 

LSS data, probing smaller scales, provide complementary information. In this paper we 
consider the real space galaxy power spectrum from IRAS Point Source Catalogue Redshift 
(PSCz) survey [19] which probes scales in the range k = 0.01 — 300 h Mpc~^ [9, 17]. In 
order to avoid problems with the interpretation of non-linear effects, we use data only at 
scales k < 0.3 h Mpc~^. 

Further, we consider the information on the linear matter spectrum at redshift ^ ~ 3 
which can be inferred from the Ly-a forest in quasar spectra. The fact that the nonlinear 
scale is smaller at higher redshift makes it possible to probe the linear power spectrum to 
smaller scales {k ~ 0.2 — 5 h Mpc~^) than are accessible to galaxy surveys at low redshift. 
In this paper we use a recent determination [10] of the linear matter power spectrum from 
the Lyman-o; forest at redshift z = 2.72, based on a large sample of Keck HIRES and Keck 
LRIS quasar spectra. To convert the measured power spectrum of Ly-a flux into the matter 
power spectrum, [10] apply a correction factor obtained from TV-body computer simulations. 
To allow for possible systematic uncertainty [20] in this correction factor, we will repeat the 
analysis excluding the data at the smallest scales. 



2.1 Data analysis 

In order to investigate how the CMB, PSCz and Ly-a data constrain the SR parameter 
space (n, r, Sinjfc), we performed a likelihood analysis of the data sets from COBE [18], 
Boomerang [15] and MAXIMA [16], together with the decorrelated linear power spectrum 
of PSCz galaxies for k <0.3 h Mpc"^ [9, 17], and the Ly-a data from Table 4 of ref. [10]. 
The likelihood function is 

/:ocexp(-xV2), (3) 

where 

2 2 I 2 I 2 /A\ 

X = XCMB + XPSCz +XLy-a- (4) 

For the CMB data 
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while for PSCz and Ly-a data 



XpSCz,Ly-a 



{Pk,i{G) — Pk,i)'^ 



(6) 



the sum being taken over published values of band-powers Ci^i and Pk^i- The quantity 9 is 
a vector of cosmological parameters, taken here to be 



The parameters are: the matter density J7„,; the baryon density the Hubble parameter 
Hq\ the optical depth r to reionization; the normalizations QcMBj Qpsczj and QLy-a of 
the CMB, LSS, and Ly-a power spectra; and the inflationary parameters n^r^d[-ak- We 
have assumed that the universe is flat, f^A = 1 ~ as predicted by standard inflationary 
models. For all the figures we marginalize over all other parameters. We consider only two 
values for the baryon density, Q,^,h? = 0.019 and 0.030, however, as we will see, the results 
are very similar in those two cases, suggesting that the results will be similar if allowing 
^l)h? as a free parameter. 

The CMB, LSS, and Ly-a data are all subject to uncertainties in their overall normal- 
izations. The CMB groups quote estimated calibration errors for their experiments, which 
we account for by allowing the data points to shift up or down, by 10% for Boomerang [15], 
and by 4% for MAXIMA [16]. Because of uncertainty in the linear galaxy-to-mass bias for 
PSCz galaxies [17], we conservatively treat the normalization Qpscz as an unconstrained 
parameter. Ref. [10] quotes uncertainties in the overall normalization of the matter power 
spectrum inferred from the Ly-a forest, but these uncertainties are based on simulations 
with d\^k = 0; so to avoid possible bias we again conservatively treat the normalization 
Qhya as an unconstrained parameter. 

For simplicity all data points have been treated as uncorrelated in the likelihood func- 
tions, eqs. (5,6). For the CMB data, correlations between estimates C/^j of angular power 
at different harmonics I are induced by finite sky coverage, but in practice the CMB teams 
quote band-powers at sufficiently well-separated bands of / that the correlations are prob- 
ably small. For the PSCz data, the published band-powers are explicitly decorrelated. For 
the Ly-a data, the covariances between estimates of the flux power spectrum are small, 
according to Fig. 12 of [10], and this may translate into small statistical covariances in the 
inferred matter power spectrum. As mentioned earlier, it is not completely clear how good 
this translation from flux power to matter power is, and we leave this question for future 
investigation. 

We have chosen the pivot scale in Eq. (1) as Uq = 0.05 Mpc^^. This choice is made for 
convenience, since = 0.05 Mpc^^ is the scale at which wave-numbers are normalised in 
the CMBFAST code. Our results are independent of the value of Uq. 



Q — {^mi ^Aj ^bi Hq, t , 



Q,n,r,ainA:} • 



(7) 
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Figure 2: The 1 and 2 a allowed regions for the two slow roll parameters n and di^ik- The left 
panels assume a BBN prior on Jlj/i^ = 0.019, whereas the right panels are for Jlj/i^ = 0.030, 
the value which best fits the CMB data. The top row is for CMB data alone, the middle 
row is for Lyman-a data alone, and the bottom row is for the combined analysis. 
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2.2 Results 



The analysis of the constraints from CMB data was presented in detail in ref. [7] (see 
also [21]). Here we extend the analysis in [7] by including the reionization optical depth, r, 
as a free parameter. This is potentially important since there is a well-known degeneracy 
between r and the scalar spectral index n. However, the analysis turns out to prefer models 
with T a; 0, so including r leaves the main conclusions of [7] essentially unchanged, as 
can be seen in the upper panels of Figs. 2 and 3. We recall here the main conclusions 
from [7]: (1) if we allow the primordial power spectrum to bend, dink 7^ 0? then CMB data 
do not constrain the tensor to scalar perturbation ratio r; (2) if we assume a BBN prior, 
Qih? = 0.019 [24], then CMB data favour a negative bend, d[-ak < 0, corresponding to a 
bump-like feature centered at scales k ~ 0.004 Mpc^^. 

Let us now discuss the information on the inflationary parameters which emerges from 
PSCz and Ly-a data. The first observation is that PSCz data do not provide relevant 
constraints in the {n^di^k) plane. This is easily understood, because the PSCz data at 
large scales (say k < 0.03 /iMpc^^) have large errors and therefore play little role in the 
evaluation. This implies that PSCz data effectively span only one decade in k. Considering 
that the overall normalization of the data is taken as a free parameter, it is evident that 
one cannot obtain strong constraints from such a small range of scales. 

The situation is quite different with Ly-a data, which have small error bars and span 
almost two decades in k. This is shown in Fig. 2, where we present the allowed regions 
corresponding to 1 and 2 a (we define the Icr region as < Xmin + 2-3 and the 2cr region 
as < Xmin + 6-17) for both Uth'^ = 0.019 as suggested by BBN (left column), and 
Clljh^ = 0.03 as suggested by CMB (right column). The top graphs are from CMB data 
alone, the middle graphs are obtained using Ly-a data alone, and the bottom graphs are 
from the combined analysis. 

It is straightforward to understand how Ly-a data select the allowed regions shown in 
the middle panels of Fig. 2. At the small scales probed by Ly-a data, the theoretical linear 
matter power spectrum P{k) is roughly proportional to 

P{k) (x\n^{ak/rh)k'''-^ (8) 

where a = 2.205, F = flmh and k is expressed in Mpc~^ [25]. The parameter n' = 
dln'P{ka)/dlnk + 1 is the effective spectral index of primordial density perturbation at the 
scale kot ~ 3.7{Qmh'^ /O-^)^^"^ Mpc^^ representative of Ly-a data ^. On the other hand the 
Ly-a data, as discussed in [10], are well fitted by a power law, P{k) ~ A,'', with spectral 
index v = —2.47 it 0.06. This means that 

If we consider that ^ T/i < 1, the previous expression can be directy translated into 
a bound for the effective spectral index of primordial density perturbations. We obtain 
n' ~ 0.75 — 1.1, which roughly corresponds to what is shown in the middle panels of 
Fig. 2 for dink = 0. It is also easy to understand the observed correlation between n and 
dink- The effective spectral index n' can be expressed, as a function of n and di^k, by 
n' = n + ln{ka/ ko)dink- It is thus evident that the allowed region in the plane (n, di^k) lies 
roughly along lines of constant n'. 



^ The observational units for wavenumbers are (km s ^) ^ . The conversion to Mpc ^ is model-dependent 
since it requires the evaluation of the Hubble constant H{z) at redshift z = 2.72. 
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It is important to note that the regions constrained by Ly-a data are "orthogonal" to the 
regions constrained by CMB data. This means that the combined analysis (CMB+Ly-a) 
gives much stronger bounds than any of the two alone. Moreover, since CMB data provide 
a bound on fim^^i the degeneracy between n' and Th in Eq. (9) is removed, and therefore 
the Ly-a observation of P{k) oc ^-2-47±0.06 ^^^j^ directly translated into a constraint on 
n'. This is clearly shown in the lower panels of Fig. 2, from which one obtains the following 
conclusions: (1) the combined analysis strongly indicates that the bend of the primordial 
power spectrum is close to zero, being —0.02 < dink < 0-01 at the 2a level; (2) the spectral 
index n is constrained to be fairly close to n = 0.9. One notes that the limits obtained do 
not crucially depend on the assumed values of Qi,, indicating that the inclusion of Lyman-a 
data avoids the problem with CMB data alone, that different Q/, extend the allowed range 
of the spectral index beyond n « 1 [21]. However, one should note that, due essentially 
to CMB data, the goodness of the fit is quite sensitive to the assumed value of il^, the 
being substantially smaller for high Slf,. 

n^ii^ = 0.019 n^h^ = 0.030 




0.60.81.01.21.4 0.60.81.01.21.4 



n n 

Figure 3: The 1 and 2 a allowed regions in the (n, r) plane. The left panels assume a BBN 
prior on Slft/i^ = 0.019, whereas the right panels are for ri^/t^ = 0.030, the value which best 
fits the CMB data. The top panels are for CMB data alone, while the bottom panels are 
for the combined analysis. Hybrid models are to the right of the full line. 

In Fig. 3 we present the constraints obtained for the remaining SR parameter r. Specif- 
ically, we show the 1 and 2 a allowed regions in the (n, r) plane, for both low and high Q),. 
The top panels are from CMB data alone, while the bottom ones are from the combined 
analysis. It is clear that in the combined analysis r is constrained to be smaller than about 
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0.3 at 2 a. This is substantially different from the analysis of CMB data alone, where no 
constraints on r could be obtained, since a strong bend, e.g. di^k = —0.2, could allow the 
tensor component to be big [7]. It is worth noting, that Ly-a data don't probe r directly; 
they fix n to be less than 1, and dink close to zero. In this way strong constraint on r 
can be obtained from CMB data. This result is extremely important when we compare to 
theoretical models, since different classes of models predict different relationships between 
the scalar spectral index n and the tensor to scalar ratio r (see Fig. 1). The combined 
analysis, indicating n « 0.9 and r < 0.3, favours small field models and seems to exclude 
hybrid models. This is clear from the lower panels of Fig. 3, where the big n, big r (often 
referred to as hybrid models, to the right of the full line), are excluded at 2a. ^ 

An important caveat to this analysis concerns possible systematic uncertainties in the 
inference of a linear matter power spectrum from the Ly-a data. What [10] measure directly 
from observations is the power spectrum of transmitted Ly-a flux. The conversion to a linear 
matter power spectrum involves a fairly large, scale-dependent correction which [10] extract 
from collisionless computer simulations, with the Ly-a optical depth taken proportional to 
a certain power of the dark matter density [13]. It has been suggested [20] that the relation 
between baryonic and dark matter densities could introduce significant uncertainty in the 
flux-to-mass correction at small scales. Specifically, pressure effects cause the baryonic 
density Si, to be smoothed compared to the dark matter density S^. an effect that can be 
parametrized as SbCk) = exp[— (A;/A;j)^](5rf(k), with comoving filter scale ~ 35/iMpc^^ 
[14]. The procedure considered by [20] is to treat the filter scale A;j as a free parameter 
constrained only by the shape of the power spectrum of Ly-a flux. However, [10] argue that 
treating kf as a free parameter is overly pessimistic, and that if kf takes values suggested 
by hydrodynamic simulations, then the effect on the power spectrum is minor. 

To allow for the possibility that the systematic errors are underestimated at small scales, 
we repeated the analysis neglecting the last 3 data-points from the Ly-a data [10], both for 
dink = and for di^k free. The upper panels in Fig. 4 show the case where dink = 0- Here, 
the results are essentially unchanged by the removal of the data points. However, the lower 
panels show the full case. Here, the tight constraint on r disappears completely, and the 
constraint on n is significantly weakened. The reason is that the small scale data points 
are the most important for constraining dink- When these points are removed, dink is not 
nearly as tightly constrained as before, and as seen for the case where only CMBR data is 
used, a large r can be compensated by a negative bend of the spectrum. 

Thus, the very tight constraint derived above depends on the correctness of the small 
scale Ly-a data. Therefore, it is highly desirable that a better understanding of the possible 
systematic errors in determining matter power spectra from Ly-a forest observations is 
developed. 

Another worry is that the error ellipses for CMB and Ly-a are only marginally over- 
lapping. This might indicate that the two data sets are mutually inconsistent, and that 
other physical effects should be taken into account. For the Ly-a data one might think of 
massive neutrinos or warm dark matter, where both would suppress power on small scales, 
potentially allowing bigger n or bigger positive dink- 

^We have used k = 5 (sec Eq. 2), and for larger k the hybrid models move to bigger r. It is worth 
pointing out, that the classification of hybrid models [11] as in Fig. 1, and as the models to the right of 
the full line in Fig. 3 is oversimplified. In this paper we follow ref. [6], and by "hybrid models" we refer to 
potentials for which e < rj (see appendix for details). When considering more general potentials, or F-term 
hybrid inflation [12], more complicated behaviour results, and a simplifled classiflcation is impossible. 
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Figure 4: The 1 and 2 a allowed regions in the (n, r) plane. These results are obtained from 
the combined analysis, neglecting the last 3 data-points from [10]. The left panels assume a 
BBN prior on r^fe/i^ = 0.019, whereas the right panels are for rij/i^ = 0.030, the value which 
best fits the CMB data. The top panels assume d\^k = 0; whereas the lower panels have 
dink as a free parameter. Hybrid models are to the right of the full line. 



3 Conclusions 

We have considered data from both CMB and LSS to place constraints on the parameters of 
single field slow roll inflationary models. We have found that, by combining CMB data with 
the power spectrum inferred from Lyman— a forest, one obtains strong constraints on the 
inflationary parameters. We obtain 0.8 <n< 0.95, r < 0.3 and —0.02 < dink ^ 0-01 at 2a 
level. In the language of SR this means e < 0.03 and rj < 0.6, with the best fit model being 
small fields. This still leaves a large part of the SR parameter space open, but seems to 
exclude big n, big r models, often referred to as hybrid models. These constraints are much 
stronger than those from CMB data alone, and arise because the error-ellipse from Ly-a 
data is almost perpendicular to the error-ellipse from CMB data. Let us repeat, that we 
in this analysis for simplicity have assumed the absence of appreciable covariances between 
the reconstructed Ly-a mass powers. 
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A Notation 

We use the notation: rj = ae, where 

e = — , r] = M^ — and P = ^ 

see [7, 21] for details. The notation with r] = M'^V" /V, used e.g. in [2], simply corresponds 
to the substitution a — )■ a + 1. The 3 classes of SR models are small fields {a < —1), large 
fields (—1 < a < 1), and hybrid models (1 < a). One finds 

r = 2Ke, ^infe = -2^^ + 8e^(2Q!-l) and n - 1 = 2e(Q! - 2) . (10) 
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